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Про інформативність показників 
інтервального аналізу сигналів 1 зображень 


У статті розглянуті питання побудови інтелектуальних пристроїв і систем моніторингу стану об'єктів 
з урахуванням компактного кодування сигналів та зображень, зроблений аналіз оптимізації процесів 
їх введення та оперативної обробки. Описані теоретичні основи отримання інтервальних характеристик 
шумів та їх використання для визначення множинних оцінок інформативних параметрів введених 
сигналів і зображень. 


Вступ 


Сучасний розвиток інтелектуальних сенсорів та відеосенсорів, пристроїв, систем 1 
мереж тривалого моніторингу станів об'єктів різноманітної природи і призначення 
грунтується на розробці ефективних методів і алгоритмів оперативної обробки сигналів, 
зображень та відеоданих в місцях виникнення інформаційних потоків, тобто безпосе- 
редньо на об'єктах. При цьому важлива мінімізація первинних потоків даних за рахунок 
виявлення та компактного кодування інформативних 1 достовірних вхідних даних на етапі 
введення (відбору) моніторингової інформації засобами об'єктних пристроїв і систем. 
Слід зазначити, що об'єктні пристрої 1 системи, як правило, працюють в умовах багатьох 
обмежень, включаючи обмеження на продуктивність обчислювальних ресурсів, об'єм 
оперативної пам'яті, час обробки даних при наявності вхідних шумів 1 промислових 
завад. Тому важливою проблемою при побудові інтелектуальних об'єктних пристроїв, 
систем і мереж є реалізація оперативної багатофункціональної обробки сигналів та 
відеосигналів |1|, |2| з використанням високопродуктивних сигнальних процесорів, 
програмуючих логічних схем, процесорних модулів та мікроконтролерів. Ключовим 
завданням багатофункціональної обробки даних на об'єктах є виявлення найбільш 
інформативних ділянок сигналів та відеосигналів, відліки яких кодуються найбільш точно, 
з підвищеною частотою дискретизації та максимальною кількістю біт. Неінформативні 
ділянки сигналів та відеосигналів, включаючи ділянки, які вражені шумами, кодуються 
менш точно 1 більш стисло. 

Основою визначення найбільш інформативних ділянок сигналів і зображень є 
оперативне визначення динамічних та статистичних показників огинаючих сигналів 
і відеосигналів |З | та показників інтервальних методів аналізу даних |4-8|, які дозво- 
ляють отримувати множинні оцінки інформативних параметрів істинного сигналу 
навіть для ділянок з високим рівнем шуму. 

Таким чином, метою статті є розробка методології отримання замкнутих та 
компактних множинних оцінок сигналу з шумом на основі інтервальних харак- 
теристик шуму та розгляд можливостей подальшого використання такого підходу 1 
його результатів до вирішення задачі стиску даних вимірювань з мінімальними 
контрольованими втратами точності. 
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Оптимізація процесів введення та оперативної обробки 
сигналів і зображень 


При введенні даних від об'єктів важливо контролювати умови введення сигналів 
і зображень, опосередковано визначати вхідне співвідношення сигнал/шум |2), |81. 
Це дозволить контролювати ступінь спотворення вхідних даних шумами і завадами. 
Відповідно масиви первинних даних доцільно супроводжувати службовою інфор- 
мацією, яка підтверджує якість введених даних та дозволяє оптимізувати обробку 
вхідних даних: чисті від шумів фрагменти та ділянки сигналів і зображень фільтруються 
та кодуються більш якісно, а дані з шумами фільтруються спрощеними методами, 
проріджуються і кодуються менш якісно і більш стисло. Слід зазначити, що сам процес 
фільтрації сигналів і зображень спотворює первинні дані, достовірність яких суттєво 
залежить від вхідного співвідношення сигнал/шум. Тому на ділянках (фрагментах) 
сигналів та зображень з шумами молодші двійкові біти є недостовірними і немає 
смислу виконувати складні і точні операції фільтрації та стиску даних з шумами. 
Мінімально допустима частота дискретизації сигналів вибирається адаптивним чином в 
широких межах, згідно з умовою /х З /от 5 2К. Х / пах» ДО Ко 28-10 -- коефіцієнт 
підвищення частоти дискретизації сигналу /зк За Котельниковим, значення якого 
сугтєво залежить від метрологічних вимог до пристроїв введення і обробки сигналу, 
Луах - Максимальна частота сигналу. На відміну від вимірювальних сигналів потік даних 
від відеосенсорів визначається їх роздільною здатністю М х М (М - кількість пікселів у 
рядку поточного кадру, М - кількість рядків у кадрі), кількістю біт, виділених для 
кодування яскравості пікселя, а також залежно від частоти отримання кадрів і виду 
зображення (кольорове, монохромне). 

Найпростішим способом адаптивного введення сигналів є їх дискретизація з 
максимальною частотою / ух 7 2. Х / хх З Наступним прорідженням відліків сигналу 


на величину К ; - (АХ Ф у, де) - 1,2,3, 4,... - поточна величина коефщієнта прорідження 
частоти /? пах - Для простоти кодування величина /) вибирається із величин | - І, 2,4,8... 
З метою зменшення первинних інформаційних потоків 1 отримання достовірних відліків 
сигналів доцільно використовувати сигма-дельта АЦІП. Оптимальним вирішенням проб- 
леми введення сигналів є тактування роботи аналого-цифрового перетворювача з 
адаптивним періодом кодування /, 2 /(Коптіп ІС/ що АХ. ф ), де Коттіп - МІНіМаЛьЬНО не- 


обхідний коефіцієнт стиску даних, |с/ш |, - відносний показник вхідного співвідношення 


сигнал/шум, ДХ. р зб4 З -Х Ф, - поточне значення приросту профільтрованого сигналу, 
визначеного в процесі реалізації оперативної фільтрації, наприклад, ковзким способом з 
мінімальним вікном усереднення / У 4-5 відліків. З метою мінімізації інформацій- 
них потоків без втрат точності кодування відліків сигналу на ділянках з високочас- 
тотними шумами і на пологих ділянках коефіцієнт прорідження частоти дискретизації 
сигналу К; вибирається максимально можливим, а кількість достовірних біт суттєвих 
відліків при компактному кодуванні 45 - 9, І, 11, 12 вибирається залежно від ступеня їх 


«зашумленості», тобто 45 - "Цс/ш|. ). У найбільш простому випадку інформація про 
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ступінь «зашумленості» ділянок сигналів характеризується різницею ДХ/" зх НоХ 0 


між поточними значеннями вхідних (з шумами) Х/" і профільтрованих Х Ф відліків 


сигналу. Шляхом порівняння величин ЛХ/" з відповідними пороговими величинами 
В,Р.,.. Р, визначаються оцінки поточного співвідношення сигнал/шум, де 5 - кіль- 
кість рівнів (станів) величини |с/ш/|,,. 


Компактне кодування відліків сигналів грунтується на оперативному обчисленні 
величин ДАХ 7, АХ і виборі величин Лат та 49 - Для реалізації на об'єктах компакт- 


ного кодування сигналів і зображень без залучення потужних процесорів ефективним є 
використання алгоритмів фільтрації-стиску даних моніторингу, які оптимізовані за 
швидкодією і точністю та грунтуються на «сигнальному» підході, який передбачає 
попереднє згладжування кривої обвідної сигналу чи фрагмента зображення (рядка, 
стовпця, відповідної поточної кількості рядків і стовпців) та пошуку суттєвих відліків 
(СВ), до яких відносяться екстремуми та точки перегину (зміни опуклості) кривої. Всі 
інші відліки, які знаходяться між сусідніми СВ, класифікуються 1 кодуються як несуттєві 
відліки. Для досягнення необхідного К.,, пуд НА ділянках обвідної з високочастотними 
складовими точками перегину можна знехтувати, оскільки вони інформативні при аналізі 
нединамічних (гладких) ділянок обвідної. З метою експрес-аналізу вибірок сигналів та 
фрагментів зображень доцільно використовувати різні показники інформативності відпо- 


відних СВ, включаючи АХ -ХОРоХОЇ, є «Аа он та інші відносні 


показники, де Х/" - амплітуда і-го СВ, 2 - показник інформативності і -го екстремуму, 


- величина, яка характеризує інформативність і -го екстремуму по від- 


ношенню до (7 - 1)-го екстремуму; - величина, яка характеризує ін- 


Ха - я 


формативність 1-го екстремуму по відношенню до (і -- 1)-го екстремуму. 


е 
Мб 


Інтервальні характеристики сигналів та їх використання 
для обробки вимірювальних даних 


Враховуючи обмеженість по величині значень реальних сигналів, при обчис- 
ленні їх статистичних характеристик доцільно враховувати параметри цих обмежень. 
В класичних моделях сигналів теорії ймовірності дана інформація не враховується. 
Наприклад, нормальний закон розподілу, який найчастіше використовується при 
обробці даних різноманітних досліджень, не передбачає такої обмеженості. В |4-8| 
запропоновані методи та алгоритми інтервального (множинного) оцінювання інфор- 
маційних параметрів масивів вимірювальних даних з використанням апріорної 
інформації, яку бажано використати. Таким чином із обмеженого сигналу випливає 
обмеженість його середніх значень на скінченних інтервалах обробки сигналів. Дані 
інтервальні обмеження теж, як апріорна інформація, можуть використовуватись для 
отримання інтервальних оцінок інформативних параметрів сигналів та масивів вимі- 
рювальних даних. Використовуючи дану методику обробки даних в задачах іден- 
тифікації 1 оцінювання параметрів сигналів на основі виміряних параметрів отримують 
множинні оцінки в просторі цих параметрів. 
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Нехай відомо, що результати вимірювань чи спостережень у вигляді послідов- 
ності «виходів» можна представити у вигляді 

У "ХХ їй» пе23,..М, (1) 

де у, - задана скалярна числова послідовність, отримана в результаті вимірювань 

чи спостережень, х, - невідома послідовність істинних значень вимірюваних чи 


спостережених процесів, а /, - невідома обмежена послідовність, що відображає 


похибки вимірювань чи спостережень на інтервалі п - 1,М (М в соп5ї » 0), 
Процедура обробки отриманих значень для обчислення оцінок істинних зна- 
чень сигналу суттєво залежить від прийнятих припущень щодо завад вимірювань чи 
спостережень. На відміну від імовірнісного підходу, розглянемо випадок, коли вико- 
ристовуються інтервальні характеристики обмежених збурень |4-7|. Найпростішими 
і фізично обгрунтованими є припущення про обмеженість самих похибок (завад) та 
швидкості їх зміни (61: 
ГЯ «б зсоп5і, А «у зсоп5бі,Мп»д0, ДР еро-/ (2) 
При цьому вважається, що /, не містить постійної складової (центрована зава- 


да), тобто є чисто шумовою складовою похибок. За наявності постійної методичної 
похибки вона відноситься до невідомих істинних значень вимірюваного процесу, а її 
виділення з них є окремою задачею, яка не буде розглядатися в даній статті. Із об- 
меженості /, випливає відсутність постійної складової в Д/, . 

Розглянемо також можливість гладження шумової складової похибок та швид- 
кості їх зміни ковзними вікнами з вибраною фіксованою шириною |9)|. Саме просте 
ковзне вікно, яке використовується для гладження - це прямокутне вікно шириною 


М, 272М31, (3) 
де М 20 - деяке ціле число. При цьому формально вводиться умова, що 
ДдермМікі2,..М, Дін хо У 12-1,2,..М. (4) 


Тоді значення згладженої послідовності // визначаються за формулою 


з із М 
даму х У Д., УМ -0,1,2,..М,, пе123,М, Му М 1/2). | (5) 
із-М 
Відповідно до |4-7| це дозволяє додатково до (2) використовувати для оціню- 
вання з врахуванням (3) інтервальні характеристики виду 


т, (М) « гора «т, (М) У М «0.1,2,..., М,, п 1,2,3,..., М, (6) 
а також аналогічні ернатній характеристики швидкості зміни завад вимірювання 
Дт, (М) « Щ ум « Дт, (М) У М « 0,1,2,..., М,, п 12,..., М - 1, (7) 
де | 
Мо КМ -2)/3), (3) 


ти( М) і т, (М) - відповідно нижня і верхня межі інтервальної функції оцінки 
арифметичного середнього значення завад, що визначається шириною інтервалу Мб, 


на якому розглядається ця оцінка, а Дт, (М) і Дт, ( /) - нижня і верхня межі анало- 


гічної інтервальної функції оцінки арифметичного середнього значення швидкості 
зміни завад. 
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Слід зауважити, що 
Д/,- М мА; -Мчі кА, -Ма2 Ж. А 


і Лема ан "Льно До Жебвя оо В р ро 
сор а. аа Уа Ро Лу ЕР п ро ЕР ду зима ее 


тобто фактично замість (6) можна написати нерівності виду 


Дт, (М) « з ня - 3. у) Є Ат, (М) У М «04,2,..., М,, пе12,..,М -1. (9) 
0 


Відзначимо, що інтервальні характеристики (6) 1 (9) є більш загальні від (2) 1 
містять їх в собі як частинний випадок при МУ-0 і т,(0)--т,(0)-6 та 


Дт, (0) - -Дт, (0) є у. Безперечно, що знання при М » 0 значень меж інтервальних 


функцій арифметичного середнього завад 1 швидкості їх змін є додатковою інформа- 
цією, яка дозволяє в більшості випадків суттєво покращити множинні оцінки відповід- 
них параметрів, ніж коли використовувати лише знання обмежень при М - 0). Але для 
отримання цієї інформації необхідні додаткові метрологічні дослідження вимірю- 
вальних пристроїв, які передбачається використовувати. А саме, слід провести серію 
вимірювань еталонних сигналів, що близькі за властивостями (насамперед, за частот- 
ним складом і динамічним діапазоном) до природних. Це дозволить виділити шумову 
складову похибок вимірювань і шляхом гладження її ковзними вікнами різної шири- 
ни виду (5) безпосередньо отримати значення меж відповідної інтервальної функції, 
що фігурують в (6). Межі інтервальної функції оцінки арифметичного середнього 
значень швидкості зміни завад отримуються перевіркою (9) при кожному вибраному 
М » 0, тобто вимогою, щоб для виділеної шумової складової 


Дт У) в во шіп а(, -Мчі ні 
Дт, (Му з ах Рама» М 2 01,2,.М,. 


Нехай істинні значення процесів залежать від деякої заданої числової послідов- 


(10) 


ності «входів» - и, (п - 1,М), причому ця залежність має вигляд 


5 
х, 7 У у хФу(и,п), пе12,.М, а1) 
Ка 


де ф, (х)- відомі функції своїх аргументів (так звані «опорні» функції), а /, (к «К - 1,5, 
5 з сопзі ) - Невідомі постійні параметри. Для зручності подальшого аналізу відомі функ- 
ції об'єднані у вектор-рядок Ф, -(Ффі(и,,п);Фо (и, п);....Фе (и, п)), а невідомі парамет- 
ри - у вектор І «( 4715 ,...318). Для визначення оцінок значень цих параметрів з до- 


статньою точністю за результатом обробки М отриманих значень у, (п - 1,М ) та отри- 


мання таким чином достовірної математичної моделі розглядуваного процесу, що дає 
можливість подальшого його прогнозування, може бути застосовано інтервальний 
(множинний) аналіз (51, (6). При цьому основою такої обробки є співвідношення 


7 У ІьхФаИн)я , пе12,М. (12) 
Ка 


Для використання обмежень, пов'язаних із швидкістю зміни завад, будуть вико- 
ристовуватись співвідношення 


5 
ду, 7 21, х АФ (и, п) А, пе12,..,М -1, (13) 
Ка 


де Дуп З Упзі 7 Уп» АФі(ип п) афи(йрапЯ1)-ФЕ(и, п). 
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Наявність завад не дає можливості в результаті обробки отримати точне значення 
параметрів, але, як випливає із обмежень (6) і (9), дозволяє отримати їх гарантовані 
оцінки у вигляді множини в просторі параметрів 3", що виділяється нерівностями, 
котрі отримуються з умови задоволення при значеннях параметрів з цих множин 
співвідношень (6) 1 (9) при всіх п і для деяких значень /, 1 Д/,, що задовольняють 
даним обмеженням. Якщо ж в результаті обробки отримаємо пусту множину, то це 
означатиме, що або гіпотеза про зв'язок істинних значень зі «входами» не виконується 
(через неправильний вибір вигляду чи кількості «опорних» функцій), або інформація 
про інтервальні характеристики завад вимірювання неточна. Такий випадок тут розгля- 
дати не будемо. 

Сформуємо процедуру обробки даних. Для цього перепишемо (12) у вигляді 


у, "Фухіеї, пе12,.,М. (4) 


Тоді, можна записати вираз для суми будь-яких рівностей (12) підряд, 
розділеної на їх кількість, починаючи з деякої (п - М)-ї до (п-- М )-ї рівності 


1 ізМ ізМ 
-УПи Фу хі) Й У М а 04,2,..М,, па М-1,М я 2,..М-М. 
Мом Мо У 
Звідси на основі (5) отримаємо систему лінійних нерівностей 
1 ізМ 
т, (М) У М ди ФІ 1) «т, (У) 
із- М 


У М «0,1,2....М,, па М41,М -2,...М - М, 


(15) 


що визначає множину значень вектора /, в просторі параметрів Е йо котрі задоволь- 
няють всі нерівності (15) при заданих межах т,(М) і т, (М) інтервальної функції 
оцінки арифметичного середнього значення завад. Систему нерівностей (15) можна 
переписати в компактному вигляді 


У(п, М) - т, (М) « Ф(п, Мух І, « У(п, М) - т, (М) 


(16) 
У М «01,2... Мі, пе М 1,М 4 2,.М М, 


де фігурують усереднені значення виміряних «виходів» і вибраних «опорних» функцій 
ізМ 


ізМ 
Упніо Фи,Муа УФ. 
М М ТУ 


Аналогічно системі (15) може бути записана на основі (13) система лінійних нерів- 
ностей 


1 
уУ(п, У) г-- 
У( Усе 


іє- 


ізМ 
дт, (М) є х У (Ду, АФ, х І) «Ат, (М) 7) 


іс-М 
М М с 0,1,2,...М,, па М 41,М -2,...М - М- 1, 
яка з врахуванням (9) може бути переписана в компактному вигляді 


Урнма 7 Унем 7 Ат, (М) « (Ф, а Фу) 02 У ун 7 Упем 7 Ат, (У) (18) 
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М М «0.1,2,...М.,,па М31,Мч2,.М -М- 1. 


Очевидно, що множина значень вектора /. в просторі параметрів Е 5 якій належать 
істинні значення невідомих параметрів, визначається об'єднаною системою лінійних не- 
рівностей (16), (13). 

Якщо відома деяка апріорна інформація про оцінку істинних значень невідомих 
параметрів у вигляді системи лінійних нерівностей відносно цих параметрів, то резуль- 
туюча множина-оцінка О; компонент вектора / отримується об'єднанням апріорної 
системи з системою лінійних нерівностей (16), (18), що визначає перетин апріорної та 
апостеріорної множин-оцінок. При цьому з неї можна отримати оцінки окремих пара- 
метрів у вигляді інтервальних зовнішніх оцінок 


і є ГІбтітад р ут що ід (13), КО з тах (1;)9 тіп (-1;), )є12,мт, (19) 
ІеО ІеО ІеО; 


що не можуть бути покращені без використання додаткової інформації, так як ці оцінки 
отримуються із рішення відповідно до (19) задач лінійного програмування. Тоді як точ- 
кові оцінки значень невідомих параметрів /; можуть вважатись серединою інтервальних 


оцінок (19). В принципі, як точкова оцінка може використовуватись будь-яка точка 
множини О;, хоча природно вибирати в цій якості деяку середню (так звану серед- 
ньозважену) точку даної множини. Отримання 1з (19) хоча би однієї пустої інтервальної 
оцінки (за рахунок суперечності виду а » нм ), означатиме, що результуюча множи- 


на-оцінка О; пуста. 

Дана методика може бути використана для досягнення оптимального значення 
коефіцієнта стиску даних при їх передачі каналами зв'язку. А саме, будемо вважати, що 
масив довжиною М даних попередньо зберігається в деякому буфері (блоці пам'яті) 
перед лінією передачі. Тоді, отримавши множинну оцінку невідомих параметрів істин- 
ного сигналу на цьому масиві, можна оцінити ступінь максимальної неточності його 
відтворення при виборі того чи іншого коефіцієнта стиску. Більше того, якщо набір 
опорних функцій і допустима неточність відтворення сигналу попередньо вже були 
вибрані, то можна попередньо розбити апріорно задану множину можливих невідомих 
параметрів на сукупність підмножин, для яких гарантується певна точність відтворення 
істинного сигналу при певному коефіцієнті стиску. Тоді треба буде тільки перевірити 
належність результуючої множини-оцінки О; до певної підмножини чи декількох під- 
множин та вибрати з них оптимальний коефіцієнт стиску. 

Такий підхід може бути поширений на випадок «матричного» сигналу на виході, 
коли зображення отримується у вигляді виходу відеосенсора на основі приладу Із 
зарядним зв'язком. Тобто замість (1) маємо 


Яка "Яра Да» УЗ МУ Бак ньо; (20) 


де Н21,021 - деякі цілі числа. Але суттєвою відмінністю тут є необхідність враху- 
вання зв'язків між значеннями істинного сигналу для різних /) 1 Кк залежно від сусідства в 


матриці вибраних елементів, які складають цілісне зображення. Також слід врахувати 
взаємозв'язок між елементами «матричної» шумової складової в (20). 


Висновки 
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Запропонована методологія обробки даних об'єктними пристроями 1 системами 
орієнтована на використання при побудові бортових систем збору, обробки 1 передачі 
інформації, в телемедичних мережах, мережах екомоніторингу та дистанційного конт- 
ролю станів промислових об'єктів, людино-машинних комплексів 1 систем. Викорис- 
тання в ній інтервальних характеристик реальних сигналів з шумами суттєво відрізняє 
її від загальноприйнятого на даний час імовірнісного підходу. Це дозволяє отримувати 
гарантовані оцінки інформативних параметрів шумів і, таким чином, підвищує досто- 
вірність розрахованих на їх основі показників масивів вибірок сигналів. За рахунок 
цього забезпечується близьке до оптимального завантаження каналів зв'язку мереж 
передачі даних. 
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М.М. Льчак, Б.М. Шевчук 

Об информативности показателей интервального анализа сигналов и изображений 

В статье рассмотрень вопросьг создания интеллектуальньгх устройств и систем мониторинга состояния 
обьектов с расчетом компактного кодирования сигналов и изображений, проведен анализ оптимизациий 
процессов их введения и оперативной обработки. Описаньт теоретические основьт получения интервальньх 
характеристик шумов и их использования для обозначения множественньх оценок информативньх 
параметров введенньх сигналов и изображений. 
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